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[i]  Water  vapor  and  methane  mixing  ratios  measured  by  the  Halogen  Occultation 
Experiment  (HALOE),  the  Atmospheric  Chemistry  Experiment  (ACE),  and  the  Michelson 
Interferometer  for  Passive  Atmospheric  Sounding  (MIPAS)  along  with  simulations 
from  the  NRL  CHEM2D  middle  atmosphere  model  are  used  to  study  the  hydrogen  budget 
in  the  equatorial  upper  stratosphere.  Multiyear  time  series  of  equatorial  upper  stratospheric 
H20  +  2*CH4  show  temporal  variations  during  periods  of  relative  long-term  stability 
in  water  vapor  and  methane  entering  the  stratosphere.  These  variations,  which  are 
anticorrelated  to  CH4,  are  quasibiennial  and  seasonal  in  nature,  and  peak  near  2  hPa  with  a 
magnitude  of  3%  of  the  H20  +  2*CH4  mixing  ratio.  We  find  that  the  ratio  of  the  changes 
in  water  vapor  to  the  changes  in  methane  is  consistently  >2  over  the  entire  equatorial 
upper  stratosphere  in  HALOE,  ACE,  and  MIPAS  data.  Ratios  of  the  H20  to  CH4  changes 
calculated  from  the  CHEM2D  model  are  similar,  but  slightly  smaller,  and  the  variations  in 
H20  +  2*CH4  are  balanced  by  nearly  equivalent  variations  in  molecular  hydrogen.  We 
use  this  relationship  to  infer  molecular  hydrogen  mixing  ratios  from  the  observations 
which  show  that  its  mixing  ratio  decreases  with  altitude  above  5  hPa.  This  net  loss  in 
molecular  hydrogen  drives  additional  water  vapor  production  so  that  multiyear  average 
profiles  of  H20  +  2*CH4  from  HALOE,  ACE,  and  MIPAS  show  an  ~0.4  ppmv  increase 
between  5  hPa  and  1.5  hPa,  so  that  total  hydrogen  is  conserved.  Collectively,  these  results 
illustrate  the  importance  of  molecular  hydrogen  in  the  equatorial  upper  stratospheric 
hydrogen  budget. 
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equatorial  upper  stratosphere,  J.  Geophys.  Res.,  115 ,  D04302,  doi:10.1029/2009JD012135. 


1.  Introduction 

[2]  The  hydrogen  budget  of  the  stratosphere  is  dominated 
by  water  vapor,  methane,  and  molecular  hydrogen  so  that 
total  long-lived  hydrogen,  Hy,  can  be  adequately  defined  as 
the  summation  of  these  three  trace  gases.  Water  vapor, 
methane,  and  molecular  hydrogen  enter  the  stratosphere 
primarily  through  the  tropical  tropopause  region,  where  the 
seasonal  variation  in  entry  water  vapor  mixing  ratios  is 
related  to  the  seasonal  variation  in  tropical  tropopause 
temperatures  [ Newell  and  Gould-Stewart,  1981;  Mote  et  ah, 
1996;  Sherwood  and  Dessler,  2000;  Holton  and  Gettelman, 
2001],  Oxidation  of  methane  and  molecular  hydrogen  create 
a  source  of  water  vapor  in  the  upper  stratosphere  and  to  a 
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lesser  extent  in  the  lower  and  middle  stratosphere  [Brasseur 
and  Solomon,  2005].  The  net  result  of  these  processes  is  to 
create  a  global  stratospheric  distribution  of  water  vapor  with 
mixing  ratios  which  increase  with  increasing  altitude  and 
latitude  [Remsberg  et  al,  1984;  Harries  et  al,  1996b]. 
Water  vapor  affects  the  stratospheric  radiation  budget 
[Forster  and  Shine,  2002],  influences  the  chemical  budget 
of  ozone  [Evans  et  al,  1998],  and  is  a  tracer  of  dynamical 
processes  due  to  its  long  photochemical  lifetime  [Hall  and 
Waugh,  1997;  Mote  et  al,  1998], 

[3]  The  chemical  process  of  methane  and  molecular 
hydrogen  oxidation  involves  a  lengthy  series  of  reactions, 
but  is  believed  to  be  reasonably  well  understood  [e.g.,  Le 
Texier  et  al.,  1988;  Rockmann  et  al.,  2004].  Briefly,  methane 
is  oxidized  by  several  chemical  reactions  to  produce  water 
vapor  and  formaldehyde  (CH20).  CH20  is  then  quickly 
destroyed  both  by  photolysis  which  yields  molecular 
hydrogen,  and  by  oxidation  to  produce  water  vapor.  Finally, 
molecular  hydrogen  can  also  be  oxidized  (in  a  similar 
manner  as  CH4)  into  water  vapor.  The  exact  amount  of 
water  vapor  produced  by  these  processes  ultimately  depends 
on  the  difference  between  the  production  of  molecular 
hydrogen  from  CH20  photolysis  and  the  loss  of  molecular 
hydrogen  through  oxidation  and  varies  somewhat  with  time 
of  season  and  location  in  the  stratosphere.  If  these  two 
effects  are  balanced,  then  molecular  hydrogen  will  maintain 
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a  constant  mixing  ratio  and  the  quantity  H20  +  2*CH4  will 
also  remain  constant.  In  this  case,  two  molecules  of  H20  are 
produced  for  each  CH4  molecule  lost,  and  the  ratio  of  the 
changes  of  H20  to  CH4  will  also  be  2.  If,  however,  the 
production  and  loss  terms  for  molecular  hydrogen  are  not 
balanced,  then  H20  +  2*CH4  will  vary,  and  the  ratio  of  H20 
to  CH4  changes  will  deviate  from  2. 

[4]  Numerous  investigations  of  the  stratospheric  hydro¬ 
gen  budget  using  H20,  CH4,  and  H2  data  have  taken  place 
shedding  detailed  information  of  the  relative  contributions 
of  the  three  gases.  In  particular,  several  studies  of  the  ratio 
of  the  changes  of  H20  to  CH4  carried  out  in  the  lower  and 
middle  stratosphere  at  various  latitudes  have  shown  values 
near  or  slightly  less  than  2  so  that  molecular  hydrogen  is  in 
a  steady  state  and  H20  +  2*CH4  is  nearly  conserved 
[Dessler  et  al.,  1994;  Engel  et  al.,  1996;  Abbas  et  al., 
1996;  Zdger  et  al,  1999;  Hurst  et  al.,  1999;  Nassar  et  al., 
2005],  Studies  with  data  from  the  Limb  Infrared  Monitor  of 
the  Stratosphere  (LIMS)  and  the  Stratospheric  and  Meso¬ 
spheric  Sounder  (SAMS)  on  Nimbus  7  and  from  HALOE 
have  revealed  significant  deviations  in  the  ratio  of  the 
changes  of  H20  to  CH4  from  2.0  in  the  equatorial  and 
midlatitude  upper  stratosphere  [Remsberg  et  al.,  1984; 
Hansen  and  Robinson,  1989;  Remsberg  et  al.,  1996]  with 
values  significantly  greater  than  2.  These  studies,  while 
motivating  further  study,  were  done  with  limited  data  or 
during  periods  of  significant  trends  in  water  vapor  which 
might  have  resulted  in  an  incorrect  calculation  of  the  ratio. 
Modeling  results  from  Le  Texier  et  al.  [1988]  show  that  the 
ratio  of  the  changes  of  H20  to  CH4  exceeds  2  in  the  upper 
stratosphere  when  the  H20  source  from  net  H2  loss  is  taken 
into  account.  Harries  et  al.  [1996a]  inferred  H2  in  the  upper 
stratosphere  and  mesosphere  from  HALOE  VI 9  data  based 
on  nonconservation  of  H20  +  2*CH4.  However,  we  will 
show  that  the  Harries  study  used  HALOE  H20  and  CH4 
measurements  which  have  biases  in  the  upper  stratosphere 
resulting  in  an  underestimate  of  H2  losses  in  this  region.  We 
also  note  a  recent  modeling  study  by  Riese  et  al.  [2006] 
which  showed  that  estimated  increases  of  tropospheric  H2 
of  ~0.55  ppmv  from  a  possible  future  hydrogen  economy 
would  lead  to  a  0.3 -0.4  ppmv  increase  in  upper  strato¬ 
spheric  H20  due  to  increased  H2  oxidation. 

[5]  This  paper  uses  observational  data  of  water  vapor  and 
methane  from  HALOE  on  the  Upper  Atmosphere  Research 
Satellite  (UARS)  [Russell  et  al.,  1993],  ACE  on  SCISAT-1 
[Bernath  et  al.,  2005],  and  MIPAS  on  the  European  Envi¬ 
ronmental  Satellite  (Envisat)  [ Fischer  et  al.,  2008]  to  study 
the  hydrogen  budget  in  the  equatorial  upper  stratosphere. 
The  data  are  analyzed  over  time  periods  when  long-term 
water  vapor  changes  in  the  upper  stratosphere  are  small 
(1996-2007)  so  that  the  hydrogen  chemistry  can  be  better 
understood.  A  large  amount  of  methane  oxidation  occurs  in 
this  region,  providing  good  water  vapor  to  methane  frac¬ 
tional  response  for  this  study.  It  is  also  a  region  of  ascent, 
and  therefore  should  be  only  minimally  affected  by  the 
Lyman-alpha  photodissociation  of  water  vapor  which 
occurs  in  the  upper  mesosphere,  and  which  could  compli¬ 
cate  a  study  of  the  conservation  of  Hy  =  H20  +  2*CH4  +  H2. 
While  the  tropical  lower  stratosphere  retains  the  signature  of 
seasonal  variations  in  water  vapor  entering  the  stratosphere, 
these  variations  are  small  in  the  upper  stratosphere  due  to 
mixing  of  air  parcels  which  occur  during  transport  from  the 


tropical  tropopause  to  the  tropical  upper  stratosphere  on 
time  scales  of  several  years.  Results  from  the  three  instru¬ 
ments  highlight  the  importance  of  net  H2  loss,  in  addition  to 
methane  oxidation,  for  the  hydrogen  budget  of  the  equato¬ 
rial  upper  stratosphere. 

2.  Instruments  and  Data  Products 

2.1.  HALOE 

[6]  The  HALOE  instrument  uses  solar  occultation  to 
measure  absorption  of  solar  energy  in  several  broadband 
and  gas  filter  spectral  bands  at  infrared  wavelengths  [Russell 
etal.,  1993].  The  calculation  of  the  fractional  response  of  solar 
radiation  makes  the  instrument  essentially  self-calibrating 
and  results  in  stable  measurements  over  long  periods  of 
time.  HALOE  measurements  are  available  from  October 
1991  to  November  2005,  and  the  third  “public”  release  data 
set  (version  19)  is  used  for  this  study. 

[7]  The  H20  and  CH4  measurements  from  HALOE  have 
been  used  extensively  in  numerous  scientific  studies  related 
to  chemistry,  dynamics,  and  modeling.  The  H20  profiles 
have  a  vertical  resolution  of  slightly  coarser  than  2  km. 
Validation  studies  by  Harries  et  al.  [1996b]  confirmed  the 
version  17  H20  profiles  compare  favorably  with  indepen¬ 
dent  measurements  and  estimated  their  precision  at  ^7% 
and  accuracy  at  ~14%  in  the  upper  stratosphere  with 
increasing  values  toward  the  top  (upper  mesosphere)  and 
bottom  (near  tropopause)  of  the  profiles.  The  version  17  H20 
errors  found  from  the  Harries  study  are  similar  to  errors  for 
the  version  19  data  [Kiev  et  al,  2000]  and,  thus,  applicable 
to  this  study.  The  CH4  measurements  have  a  vertical 
resolution  of  ~4.5  km  due  to  the  larger  retrieval  vertical 
spacing  compared  to  the  H20  retrieval.  The  version  1 7  CH4 
measurements  were  validated  by  Park  et  al.  [1996]  who 
calculated  precision  and  accuracy  values  near  the  strato- 
pause  of  ~2%  and  ~6%,  respectively.  The  fractional  errors 
increase  slightly  both  toward  the  lower  stratosphere  and 
middle  mesosphere.  Version  17  CH4  errors  are  also  similar 
to  those  for  the  version  19  data  (J.  M.  Russell  III,  personal 
communication,  2009). 

2.2.  ACE-FTS 

[8]  The  ACE-FTS  instrument  (hereafter  referred  to  sim¬ 
ply  as  ACE)  is  a  high-resolution  Fourier  Transform  Spec¬ 
trometer  which  measures  atmospheric  absorption  spectra 
using  solar  occultation  between  2.2  and  13.3  /jm  (750- 
4400  cm-1)  [Bernath  et  al.,  2005].  ACE  measurements  are 
available  from  February  2004  to  the  present,  and  the  version 
2.2  data  are  used  for  this  study. 

[9]  H20  data  from  ACE  have  been  validated  by  Carleer 
et  al.  [2008]  by  comparisons  to  H20  measurements  from 
several  spacebome  instruments.  The  authors  concluded  that 
the  measurements  provided  an  accurate  representation  of 
H20  to  a  level  better  than  5-10%  in  the  stratosphere.  The 
repeatability  of  the  measurement  relative  to  correlative 
measurements  and  its  high  precision  (<5%)  make  it  very 
suitable  to  scientific  studies  and  a  good  reference  measure¬ 
ment.  ACE  CH4  data  have  been  validated  by  De  Maziere  et 
al.  [2008],  The  accuracy  of  individual  CH4  measurements  is 
within  25%  in  the  middle  and  higher  stratosphere;  however 
this  estimate  was  based  primarily  on  a  comparison  to 
HALOE  CH4  which,  as  will  be  discussed  later,  has  low 


2  of  16 


D04302 


WROTNY  ET  AL.:  HYDROGEN  BUDGET  OF  UPPER  STRATOSPHERE 


D04302 


mixing  ratios  in  the  tropical  upper  stratosphere  compared  to 
both  ACE  and  Atmospheric  Trace  Molecule  Spectroscopy 
(ATMOS)  data.  Comparisons  of  ACE  and  MIPAS  CH4 
coincident  profiles  show  differences  within  10%  in  the 
upper  stratosphere  [Hopfner  et  al,  2007;  De  Maziere  et 
al.,  2008],  suggesting  that  the  25%  accuracy  estimate  by  De 
Maziere  is  likely  too  large.  The  precision  values  for  ACE 
H20  and  CEI4  are  statistical  fitting  errors  from  the  least 
squares  based  retrieval  process  [Boone  et  al.,  2005]  and 
have  values  of  2-3%  in  the  stratosphere.  The  vertical 
resolution  of  H20  and  CH4  profiles  is  typically  near  4  km. 

2.3.  MIPAS 

[10]  The  MIPAS  instrument  is  a  midinfrared  Fourier 
transform  emission  spectrometer  which  measures  limb 
spectra  with  high  spectral  resolution  over  the  wavelength 
range  of  4.15  to  14.6  //m.  MIPAS  measurements  have  been 
made  from  September  2002  to  the  present,  but  problems 
with  the  interferometer  drive  unit  limited  operations  during 
the  early  2004  to  December  2007  time  period  [Fischer  et 
al.,  2008].  This  study  uses  data  produced  by  a  retrieval 
processor  designed  at  the  Institute  for  Meteorology  and 
Climate  Research  (1MK)  [von  Clarmann  et  al.,  2003]  and 
the  Institute  de  Astrofisica  de  Andalucia  [Funke  et  al., 
2001],  We  used  data  with  version  numbers  V30H2011 
and  V30  H2012  for  H20,  and  V30_CH4_8  and 
V30CH411  for  CH4,  with  the  higher  version  number 
being  used  to  achieve  better  data  sampling  for  the  last  few 
months  of  the  nineteenth  month  time  series  (September  2002 
to  March  2004)  used  in  the  analysis.  The  differences  between 
these  subtypes  of  V30  retrievals  is  small  (<0.1  ppmv  for 
H20  and  <0.01  ppmv  for  CH4),  and  should  have  no  affect 
on  the  results  presented  here. 

[11]  IMK/IAA  MIPAS  H20  data  have  been  validated 
versus  numerous  instruments  by  Milz  et  al.  [2005,  2009], 
and  a  comparison  to  the  ACE-FTS  instrument  is  given  by 
Carleer  et  al.  [2008].  MIPAS  H20  was  found  to  be  ~10% 
high  versus  HALOE,  while  the  consistency  with  ACE-FTS 
data  is  on  the  order  of  5%  with  MIPAS  being  dryer  than 
ACE-FTS.  A  northern  midlatitudinal  error  budget  analysis 
by  Milz  et  al.  [2005]  found  that  the  total  H20  random  and 
systematic  errors  were  less  than  10%  in  the  stratosphere. 
Our  analysis  of  the  MIPAS  H20  error  budget  has  confinned 
that  these  errors  are  also  valid  for  tropical  data  for  altitudes 
above  35  km.  CH4  has  been  found  to  be  high  by  5-10% 
compared  to  ACE-FTS  in  the  upper  troposphere  to  middle 
stratosphere  [Hopfner  et  al,  2007]  (also  M.  Engelhardt  et 
al.,  Institut  fuer  Meterologie  und  Klimaforschung,  Univer- 
sitaet  Karlsruhe,  unpublished  results,  2007).  The  accuracy 
of  a  southern  midlatitude  CH4  profile  was  estimated  to 
range  from  10%  in  the  upper  troposphere  up  to  ~20%  in  the 
upper  stratosphere  with  precision  values  of  ~10%  through¬ 
out  the  stratosphere  [Glatthor  et  al,  2005].  Our  error 
analysis  has  shown  that  the  midlatitude  upper  stratospheric 
error  values  for  CH4  quoted  above  are  very  similar  to  their 
values  at  tropical  latitudes. 

3.  Creation  of  a  Uniform  Data  Set 

[12]  Comparison  of  H20  profiles  from  HALOE  to  a 
variety  of  coincident  H20  profiles  from  other  instruments 
[Harries  et  al.,  1996b;  Kiev  et  al.,  2000;  McHugh  et  al.. 


2005;  Milz  et  al.,  2005;  Lambert  et  al.,  2007;  Nedoluha  et 
al.,  2007]  consistently  show  that  HALOE  H20  mixing 
ratios  are  biased  dry.  In  the  upper  stratosphere,  the  dry  bias 
is  a  few  percent  near  5  hPa  steadily  increasing  up  to  5-10% 
near  the  stratopause.  The  sense  of  the  altitude-dependent 
bias  is  similar  between  comparisons  of  HALOE  profiles 
to  those  from  ATMOS,  ACE,  MIPAS,  and  AURA-MLS 
[Harries  et  al,  1996b;  Kiev  et  al.,  2000;  McHugh  et  al., 
2005;  Milz  et  al.,  2005;  Lambert  et  al.,  2007;  Nedoluha  et 
al.,  2007].  A  comparison  of  HALOE  CH4  to  ATMOS  data 
[Park  et  al.,  1996]  and  ACE  data  [McHugh  et  al.,  2005] 
shows  that  the  HALOE  CH4  is  biased  low  relative  to  the 
other  two  data  sets  in  the  middle  and  upper  stratosphere 
with  the  bias  increasing  with  decreasing  pressure.  As 
mentioned  previously,  the  ACE  H20  data  have  a  high 
accuracy  (5-10%)  in  the  upper  stratosphere  [Carleer  et 
al,  2008],  while  the  accuracy  of  ACE  CH4  is  ^10%  in  the 
upper  stratosphere.  Also,  comparisons  of  water  vapor 
between  ACE,  the  Microwave  Limb  Sounder  (MLS)  aboard 
AURA,  and  the  ground-based  Water  Vapor  Millimeter-wave 
Spectrometer  (WVMS)  instrument  near  20°N  show  excel¬ 
lent  agreement  in  the  absolute  values  of  the  water  vapor 
mixing  ratio  at  50  km  [Nedoluha  et  al.,  2009].  Therefore, 
based  on  the  good  absolute  accuracy  of  the  ACE  measure¬ 
ments  along  with  the  similarity  of  the  HALOE  H20  and 
CH4  altitude-dependent  differences  relative  to  ACE  and 
ATMOS,  we  have  used  the  ACE  data  as  the  standard 
and  have  modified  the  HALOE  profiles  to  remove  the 
relative  bias.  While  we  have  chosen  ACE  as  the  absolute 
calibration  source  for  CH4,  we  do  note  that  comparisons  of 
ACE  CH4  profiles  to  those  from  the  MIPAS  instrument 
[Hopfner  et  al.,  2007]  show  a  somewhat  similar  altitude- 
dependent  bias  to  that  of  the  ACE  -  HALOE  differences  in 
the  upper  stratosphere,  which  could  suggest  that  there  is 
some  altitude-dependent  bias  in  the  ACE  data. 

[13]  The  bias  correction  has  been  established  using  equa¬ 
torial  HALOE  and  ACE  data  during  February  and  August 
of  2004  and  2005.  Of  the  four  months  per  year  in  which 
both  HALOE  and  ACE  made  tropical  (10°S-10°N)  meas¬ 
urements,  only  February  and  August  yielded  an  adequate 
number  of  profiles  (on  average,  32  measurements  from  each 
instrument  per  month)  for  a  meaningful  comparison. 
Monthly  zonal  mean  H20  and  CH4  profiles  over  10°S  to 
10°N  were  calculated  for  each  of  the  four  months,  averaged, 
and  then  the  average  HALOE  and  ACE  profiles  were 
differenced  as  ACE  minus  HALOE  divided  by  HALOE. 
The  difference  profiles  between  HALOE  and  ACE  H20 
were  then  smoothed  with  an  8  km  boxcar  average  to  remove 
small-scale  differences  which  may  be  related  to  differences 
in  vertical  resolution.  Figure  1  shows  the  profiles  of  the 
4  month  mean  HALOE  and  ACE  H20  and  CH4  data  along 
with  the  difference  profiles.  HALOE  and  ACE  CH4  profiles 
have  less  small-scale  vertical  structure  than  the  H20  pro¬ 
files,  primarily  due  to  the  lower  resolution  of  the  CH4 
profiles,  hence  no  smoothing  is  applied  to  the  CH4  differ¬ 
ence  profile.  The  HALOE  and  ACE  H20  and  CH4  differ¬ 
ence  curves  based  on  the  equatorial  data  comparison  differ 
slightly  from  the  results  derived  from  a  HALOE  correction 
based  on  a  comparison  of  HALOE  and  ACE  data  at  middle 
and  high  latitudes  from  McHugh  et  al.  [2005]  (not  shown), 
but  both  reveal  altitude-dependent  biases  which  increase  (by 
percentage)  with  decreasing  pressure  above  ~  1 0  hPa. 
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Figure  1.  Mean  water  vapor  and  methane  profiles  averaged  over  February  and  August  of  2004  and 
2005  and  10°S-  10°N  from  HALOE  and  ACE  and  their  percentage  differences,  (left)  Mean  profiles  (solid 
line,  HALOE;  dashed  line,  ACE),  and  (right)  profiles  of  their  percentage  differences  where  the 
differences  are  calculated  as  (ACE-HALOE)/HALOE.  The  water  vapor  differences  have  been  smoothed 
to  account  for  vertical  resolution  differences  between  the  two  instruments.  Unless  otherwise  indicated,  all 
HALOE  data  shown  and  used  in  this  paper  are  corrected  by  the  fraction  shown  here. 


[14]  Comparisons  of  MIPAS  CH4  data  to  ACE  [Hopfner 
et  al.,  2007;  Engelhardt  et  al.,  unpublished  results,  2007] 
show  that  the  MIPAS  data  have  small  positive  biases  of 
5-10%  in  the  lower  and  middle  stratosphere  and  on  the 
order  of  5%  in  the  upper  stratosphere.  A  comparison  of 
MIPAS  H20  to  ACE  H20  [. Hopfner  et  al.,  2007;  Carleer  et 
al.,  2008]  shows  differences  of  ^5%  throughout  the  strato¬ 
sphere,  with  MIPAS  generally  having  lower  values  relative 
to  ACE.  In  order  to  minimize  the  possible  influence  of  these 
biases  in  the  current  analysis,  we  have  chosen  to  correct  for 
the  biases  in  the  MIPAS  H20  and  CH4  profiles  by  compar¬ 
ing  them  directly  to  the  HALOE  corrected  profiles,  which 
means  that  the  MIPAS  corrections  are  mainly  based  on 
ACE.  We  chose  to  correct  MIPAS  to  the  corrected  HALOE 
(versus  MIPAS  to  ACE)  because  of  insufficient  temporal 
overlap  of  the  MIPAS  and  ACE  data  sets.  The  months  in 
which  both  HALOE  and  MIPAS  observed  the  tropics 
(10°S-10°N)  were  identified,  and  monthly  mean  H20  and 
CH4  mixing  ratios  corresponding  to  these  months  were 
averaged  for  each  instrument  as  a  function  of  pressure. 
The  mean  H20  and  CH4  mixing  ratios  for  each  instrument 
were  differenced  to  create  H20  and  CH4  bias  curves  as  was 
done  for  the  HALOE  -  ACE  comparison  (not  shown). 

[15]  The  equatorial  bias  corrections  have  been  used  to 
modify  the  HALOE  and  MIPAS  H20  and  CH4  profiles  for 
each  data  set.  All  subsequent  analyses  (unless  otherwise 
indicated)  will  be  performed  with  the  modified  HALOE  and 
MIPAS  data  sets.  Figure  2  (top)  shows  the  monthly  zonal 
mean  time  series  of  the  corrected  HALOE,  ACE  and 
corrected  MIPAS  H20  and  two  times  CH4  data  at  2.2  hPa 
averaged  over  10°S-10°N.  Monthly  averages  are  only 
shown  here  (and  elsewhere)  if  there  are  at  least  12  individ¬ 
ual  profiles  available.  The  2*CH4  time  series  show  quasi¬ 
biennial  (QBO)  variations  which  peak  near  2.2  hPa  in  the 
equatorial  upper  stratosphere  due  to  the  large  vertical 


gradients  in  CH4  [Cordero  et  al,  1997]  and  directly  relate 
to  the  observed  QBO  variations  in  H20  through  CH4 
oxidation.  An  increase  in  H20  mixing  ratios  during  the 
early  and  middle  1990s  is  also  evident  [Nedoluha  et  al., 
2003],  The  HALOE  and  ACE  H20  and  CH4  time  series 
agree  well  for  the  overlap  time  period  (2004-2005),  though 
differences  of  <0.2  ppmv  occur.  A  detailed  inspection  of 
individual  months  over  this  period  shows  that  when 
HALOE  H20  is  high  relative  to  ACE  H20,  the  HALOE 
CH4  measurement  is  low  relative  to  ACE,  and  vice  versa. 
This  anticorrelation  of  H20  and  CH4  differences  is  expected 
for  measurements  that  are  not  sampling  exactly  the  same 
times  and  locations,  and  is  likely  a  result  of  one  instrument 
observing  air  which  has  experienced  slightly  more  (or  less) 
methane  oxidation.  The  agreement  between  HALOE  and 
MIPAS  H20  and  CH4  is  generally  good,  but  also  shows 
small  differences  in  H20  and  CH4  which  are  anticorrelated 
in  the  same  sense  as  seen  with  the  HALOE  -  ACE 
comparison.  Nevertheless,  the  good  overall  agreement  in 
the  temporal  variations  between  the  three  H20  and  CH4 
time  series  at  2.2  hPa  and  other  pressure  levels  used  in  this 
analysis  (1.5,  2.2,  3.2,  4.6,  and  6.8  hPa;  corresponding 
approximately  to  44.1,  41.2,  38.3,  35.7,  and  33.0  km, 
respectively)  provides  confidence  in  the  corrections  to 
HALOE  and  MIPAS. 

[16]  While  the  H20  +  2*CH4  time  series  in  Figure  2 
(bottom)  clearly  shows  smaller  variations  than  either  the 
H20  or  the  CH4  time  series,  there  remain  significant 
variations.  Some  of  the  largest  variations  seem  to  be 
primarily  quasibiennial  in  nature  and  are  correlated  with 
the  H20  and  CH4  variations  in  time  with  H20  +  2*CH4 
mixing  ratios  low  when  CH4  is  high.  The  nonconservation 
of  the  HALOE  upper  stratospheric  H20  +  2*CH4  has  also 
been  noted  by  Geller  et  al.  [2002].  Clearly,  H20  +  2*CH4  is 
not  conserved  over  this  time  period  in  the  three  data  sets, 
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Figure  2.  Monthly  time  series  of  corrected  HALOE  (black  line),  ACE  (blue  line),  and  corrected  MIPAS 
(red  line)  water  vapor  and  (top)  two  times  methane  and  (bottom)  H2O  +  2*CH4  at  2.2  hPa  averaged  over 
10°S—  10°N.  The  H20  time  series  were  shifted  to  fit  the  plot.  Data  markers  on  the  top  of  each  plot  show 
the  times  of  monthly  observations  (at  least  12  profiles  required  per  month)  for  each  instrument  (same 
color  connotations). 


with  the  standard  deviation  of  the  mean  about  ±0.2  ppmv,  or 
3%  of  the  total  H20  +  2*CH4  mixing  ratio. 

[17]  We  further  illustrate  the  variations  in  H20,  CH4,  and 
H20  +  2*CE14  by  calculating  the  amplitudes  of  the  QBO, 
annual,  and  semiannual  cycles  as  function  of  pressure 
determined  by  a  least  squares  regression  fit  to  the  HALOE 
data  (Figure  3).  The  regression  fit  is  done  with  data  from 
1996  to  2005  over  10°S  to  10°N  and  included  constant, 
QBO,  annual,  and  semiannual  cycle  terms.  Data  from  1991 
to  1995  was  excluded  from  the  fit  to  avoid  a  possible 
influence  from  the  early  1990s  increase  in  H20.  The  QBO 
is  fit  with  two  empirical  orthogonal  functions  (EOF)  which 
were  determined  from  stratospheric  zonal  winds  above 


Singapore  and  provided  by  S.  Frith  (personal  communica¬ 
tion,  2008). 

[is]  In  the  lower  and  middle  stratosphere  (below  ~5  hPa), 
the  variations  are  primarily  due  to  the  seasonal  cycle  of 
tropical  tropopause  entry  mixing  ratios  of  H20  which 
propagate  upward  in  the  lower  stratosphere.  QBO  vertical 
transport  anomalies  vary  the  speed  at  which  the  seasonal  tape 
recorder  signal  ascends,  and  therefore  can  produce  a  QBO 
signal  in  this  region  even  in  the  absence  of  a  QBO  variation 
in  water  vapor  entering  the  stratosphere  [McCormack  and 
Siskind,  2002].  CH4  amplitudes  are  significantly  smaller 
than  H20  amplitudes  in  the  lower  stratosphere  and  thus 
contribute  little  to  the  H20  +  2*CH4  amplitudes.  In  the 
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Figure  3.  Profiles  of  annual  (solid  black  line),  semiannual  (dashed  red  line),  and  QBO  (dash-dotted  blue 
line)  cycle  amplitudes  from  a  least  squares  fit  to  (a)  HALOE  H20,  (b)  CH4,  and  (c)  H20  +  2*CH4  time 
series  from  1996  to  2005  averaged  over  10°S-10°N. 
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Table  1.  Single-Level  (3  Values  From  HALOE  Occultation 
Measurements  Over  1996-2005  and  10°S-10°N 


Pressure 

(hPa) 

Case  a: 
Reported 
Errors 

Case  b: 
2*H20/ 
0.5*CH4 

Case  c: 
0.5*H2O/ 
2.0*CH4 

Correlation 

Coefficient 

1.5 

4.03 

2.67 

5.91 

-0.65 

2.2 

3.09 

2.48 

3.72 

-0.80 

3.2 

2.78 

2.38 

3.08 

-0.86 

4.6 

2.81 

2.34 

3.03 

-0.86 

upper  stratosphere,  however,  there  are  large  seasonal  cycles 
in  H20  and  CH4  which  are  caused  by  seasonal  changes  in 
transport. 

[19]  The  tape  recorder  signal  and  all  variations  related  to 
it  decrease  with  increasing  altitude  up  to  ~10  hPa,  but  then 
the  variations  begin  to  increase  again.  Between  1.5  and 
10  hPa,  the  variations  in  H20  are  anticorrelated  with 
variations  in  CH4,  since  parcels  which  have  experienced 
more  CH4  oxidation  will  necessarily  have  more  H20.  Given 
complete  CH4  oxidation,  one  might  therefore  expect  that  the 
variation  in  H20  +  2*CH4  would  continue  to  decrease  with 
increasing  altitude  in  the  upper  stratosphere.  However,  as  is 
shown  in  Figure  3c,  this  is  not  the  case. 


4.  Correlation  Analysis 

[20]  We  examine  the  H20/CH4  relationship  by  studying 
the  correlation  of  H20  to  CH4.  We  then  define  the  negative 
of  the  slope  of  the  changes  of  H20  to  the  changes  of  CH4  as 
/ 3  [Le  Texier  et  al.,  1988], 

A[H20]  =  -0*A[CH4\.  (1) 

We  calculate  (3  using  regression  fits  both  for  H20  versus 
CH4  changes  at  one  pressure  level  (single-level  3)  and  for 
H20  versus  CH4  changes  over  a  range  of  pressures 
(multilevel  3).  The  single-level  3  values  are  particularly 
sensitive  to  assumptions  about  measurement  precision, 
while  the  multilevel  3  values  are  most  sensitive  to  level- 
dependent  systematic  errors  in  the  measurements.  Unless 
chemistry  changes  quickly  with  altitude,  the  methods 
should  give  similar  results. 

[21]  For  calculations  of  3  on  one  pressure  surface,  meas¬ 
urements  of  a  sufficient  time  period  are  necessary  so  that  the 
CH4  changes  are  large  enough  relative  to  the  measurement 
precision  error  in  order  to  calculate  a  linear  least  squares  fit 
line  to  the  H20/CH4  relationship  that  is  not  simply  depen¬ 
dent  on  the  errors  in  CH4.  Since  the  transport  time  scales 
that  cause  the  largest  variations  in  CH4  are  related  to 
seasonal  and  QBO  variations,  at  least  1-2  years  of  data  is 
suitable,  though  shorter  time  periods  are  possible  if  the 
measurement  precision  error  is  small  compared  to  the  size 
of  the  CH4  variations.  This  time  requirement  is  not  neces¬ 
sary  for  multilevel  3  calculations  since  large  CH4  changes 
occur  over  range  of  pressures.  It  is  also  important  that  the  3 
calculation  take  place  over  a  time  period  when  the  H20/CH4 
relationship  is  not  greatly  affected  by  nonconservative 
processes  to  either  H20  or  CH4,  for  example  a  long-term 
trend  in  stratospheric  entry  mixing  ratios  of  H20  or  CH4 
that  propagates  to  the  upper  stratosphere.  There  was  no  clear 
trend  in  H20  and  CH4  in  the  upper  stratosphere  between 


1996  and  2002  [Nedoluha  et  al.,  2003].  During  2001,  an 
influx  of  very  dry  air  was  observed  in  the  equatorial  lower 
stratosphere  [Randel  et  al.,  2004];  however,  the  change  in 
water  vapor  in  the  upper  stratosphere  (at  most  a  few  tenths 
of  a  ppmv)  appears  to  be  much  smaller  than  the  early  1990s 
increase  of  nearly  1  ppmv.  We  therefore  included  upper 
stratospheric  HALOE  data  from  1996  to  2005  in  the  3 
analysis.  Additionally,  we  separately  analyze  the  time 
period  2004-2005  for  a  direct  comparison  to  the  ACE  data. 

4.1.  Single-Level  Correlation  Analysis 

[22]  For  the  regression  at  one  pressure  level,  we  calculate 
the  mean  H20  and  CH4  mixing  ratios  for  all  measurements 
with  observation  latitudes  between  10°S  and  10°N  and  over 
1996  to  2005.  A  least  squares  regression  to  determine  a  best 
straight-line  fit  to  the  H20  to  CH4  changes  is  performed 
using  the  “fitexy.pro”  routine  in  the  Interactive  Data 
Language  (IDL)  \  Press  et  al,  1992],  and  the  negative  of  the 
slope  of  the  best  fit  line  is  taken  as  the  3  value  (equation  ( 1 )). 
The  fit  is  sensitive  to  the  random  uncertainties  in  both  the 
H20  and  CH4  measurements  which  for  H20  profiles  are 
taken  from  Harries  et  al.  [1996b]  and  for  CH4  profiles  from 
Park  et  al.  [1996],  Given  that  random  uncertainties  are 
inherently  difficult  to  quantify,  we  calculate  3  for  three 
different  estimates  of  random  error,  one  for  the  reported 
uncertainties  (case  a),  one  with  the  CH4  uncertainties  halved 
and  the  H20  uncertainties  doubled  (subsequently  referred  to 
as  a  4: 1  H20/CH4  random  error  ratio)  (case  b),  and  one  with 
the  H20  uncertainties  halved  and  the  CH4  uncertainties 
doubled  (a  1:4  H20/CH4  random  error  ratio)  (case  c).  We 
note  that  3  calculations  done  by  doubling  both  the  H20  and 
CH4  random  errors  (or  for  other  cases  where  the  H20/CH4 
error  ratio  reduces  to  1:1)  would  yield  the  same  3  value  as 
for  the  1:1  H20/CH4  error  ratio  case. 

[23]  Values  of  3  were  calculated  on  four  pressure  surfaces 
between  1.5  and  4.6  hPa,  and  for  the  three  cases  of  error 
ratios  discussed  above.  These  pressure  levels  were  chosen  to 
sample  the  region  of  the  tropical  upper  stratosphere  where 
the  H20  +  2*CH4  variations  increase  with  decreasing 
pressure  above  5  hPa  (Figure  3c).  For  pressures  greater 
than  ~5  hPa,  the  variations  in  methane  oxidation  (and, 
hence,  CH4  mixing  ratios)  are  small  and  there  may  remain 
some  small  seasonal  variations  in  H20  mixing  ratios  which 
originate  at  the  tropical  tropopause,  hence  calculations  of  3 
are  difficult  to  interpret.  3  calculations  are  also  not  possible 
above  1.5  hPa  because  most  of  the  CH4  rich  air  from  the 
tropical  lower  stratosphere  has  oxidized  by  the  time  it 
reaches  the  stratopause.  The  results  of  the  3  calculations 
are  summarized  in  Table  1.  The  correlation  coefficients  of 
the  H20  to  CH4  changes  for  each  pressure  are  also  shown  in 
Table  1  and  in  the  remaining  tables  in  this  article. 

[24]  In  Figure  4a,  we  show  an  example  of  the  regression 
plot  for  calculating  3  values  at  1.5  and  2.2  hPa.  The  data 
points  corresponding  to  individual  occultation  measure¬ 
ments  at  the  1.5  hPa  pressure  level  are  shown  in  blue,  while 
the  points  at  2.2  hPa  are  shown  in  red.  The  linear  least 
squares  regression  fit  lines  shown  here  (same  color  con¬ 
notations)  are  calculated  by  assuming  half  the  CH4  random 
error  and  twice  the  H20  random  error  at  this  level.  The 
actual  error  values  used  for  the  regressions  were  1.3%  for 
CH4  and  15.7%  for  H20  at  1.5  hPa  and  1.4%  for  CH4  and 
15.3%  for  H20  at  2.2  hPa.  The  3  value  was  2.67  ±  0.02  at 
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Figure  4.  Water  vapor  versus  methane  at  (a)  1.5  and  2.2  hPa  measured  by  HALOE  from  1996  to  2005 
and  (b)  1.5  and  2.2  hPa  measured  by  ACE  from  2004  to  2007.  The  data  points  and  linear  least  squares  fit 
line  to  the  data  are  shown  in  blue  for  1.5  hPa  and  in  red  for  2.2  hPa.  The  linear  least  squares  fit  line 
through  the  cumulative  distribution  of  points  at  1.5  and  2.2  hPa  is  shown  in  black.  All  data  shown  are  for 
the  latitude  range  of  10°S-  10°N.  The  slopes  of  the  linear  least  squares  fit  lines  are  shown  in  the  top  right- 
hand  comer  of  each  plot.  The  dashed  lines  have  slopes  of  —2.0. 


1.5  hPa  and  2.48  ±  0.01  at  2.2  hPa,  where  the  uncertainty  in 
/ 3  is  due  to  the  linear  fitting  error.  Given  the  large  number  of 
points,  the  fitting  error  represents  an  insignificant  compo¬ 
nent  in  the  total  uncertainty  in  (3.  For  comparison  to  the 
linear  least  squares  fit  lines,  a  dashed  line  with  a  slope  of 
—2.0  representing  the  conservation  of  H20  +  2*CH4  and 
assumed  steady  state  EE  mixing  ratio  is  also  shown  in 
Figure  4a. 

[25]  If  the  rate  of  EE  production  via  CH4  oxidation  is 
balanced  by  the  rate  of  EI2  oxidation  to  H20,  then  two 
molecules  of  EEO  will  be  produced  for  each  CH4  molecule 
lost.  In  terms  of  a  (3  value,  [3  would  then  be  equal  to  exactly 
two.  From  the  HALOE  analysis,  the  f3  values  are  all  greater 
than  2.0  (taking  into  account  the  small  fitting  error)  for  all 
pressure  levels  and  H20/CH4  random  error  ratio  cases, 
implying  that  more  than  two  molecules  of  H20  are  pro¬ 
duced  for  each  CH4  molecule  lost.  Since  the  only  other 
reservoir  of  hydrogen  is  H2,  this  suggests  that  there  is  a  net 
loss  of  H2  in  the  equatorial  upper  stratosphere  which,  in 
addition  of  CH4  oxidation,  drives  the  production  of  H20.  It 
is  also  noted  that  an  altitude  dependence  is  suggested  in  the 
HALOE  3  values,  with  higher  values  at  smaller  pressures. 
This  altitude  dependence  may  reflect  variations  in  the 
balance  between  H2  production  and  loss  as  a  function  of 
altitude. 

[26]  The  single-level  /3  values  from  Table  1  show  large 
variability  for  the  three  different  random  error  ratio  cases. 
These  3  values  are  calculated  from  individual  occupation 
measurements  of  H20  which  are  often  noisy  in  the  upper 
stratosphere.  The  calculation  of  /3  from  noisy  data  will  result 
in  a  large  variation  in  values,  depending  on  the  random  error 
estimate.  We  therefore  also  calculate  3  from  month-to- 
month  changes  in  H20  and  CH4  instead  of  from  the 
occupation  measurements  for  the  three  random  error  cases 
described  above  (Table  2).  This  averaging  will  have  the 
effect  of  removing  most  of  the  random  variability  in  the  data 
while  still  retaining  the  month-to-month  changes.  We  use 
the  same  random  errors  in  H20  and  CH4  as  for  the 
occupation  data  analysis  but  now  reduce  them  by  the  square 
root  of  the  average  number  of  HALOE  profiles  per  month. 
The  (3  values  show  much  smaller  variations  for  the  three 


H20/CH4  random  error  ratio  cases  which  results  from  an 
averaging  of  the  noisy  individual  measurements.  We  note 
that  the  linear  fitting  errors  using  monthly  means  are  similar 
to  those  found  using  the  occupation  data  (0.01-0.02  at 
1.5  hPa).  This  analysis  emphasizes  the  potential  effect  of 
measurement  noise  on  determining  the  slope  from  a  linear 
least  squares  fit  line.  Also,  for  a  given  H20/CH4  random 
error  ratio,  the  calculated  3  value  should  be  nearly  identical 
using  either  individual  measurements  or  for  data  averaged 
over  time  (e.g.,  monthly).  A  comparison  of  the  [3  values 
calculated  from  the  occupation  data  (Table  1)  and  the 
monthly  data  (Table  2)  shows  that  the  best  agreement  is 
achieved  by  assuming  half  the  CH4  random  error  and  twice 
the  H20  random  error  (case  b).  These  results  suggest  that 
the  ratio  of  H20  to  CH4  random  errors  may  be  larger  than 
have  been  previously  reported  [ Harries  et  al.,  1996b;  Park 
et  al.,  1996]. 

[27]  Correlation  plots  of  ACE  H20  and  CH4  are  made  in 
a  similar  manner  to  the  HALOE  analysis  except  the  data 
period  is  from  2004  to  2007.  Due  to  the  inclination  of  the 
ACE  orbit,  the  data  sampling  in  the  tropics  is  sparser  than 
that  of  HALOE.  Figure  4b  shows  the  correlations  of  the 
ACE  H20  and  CH4  at  1.5  and  2.2  hPa  where  a  H20/CH4 
random  error  ratio  of  1 : 1  has  been  used  for  the  regression 
fits.  The  calculated  (3  was  2.89  ±  0.03  at  1.5  hPa  and  2.61  ± 
0.10  at  2.2  hPa,  statistically  greater  than  2.0  and  persistent 
for  the  later  time  period  compared  with  the  HALOE 
analysis. 

[28]  The  (3  values  were  calculated  at  the  pressure  levels 
from  1  to  4.6  hPa  (Table  3).  We  only  show  results  where  the 


Table  2.  Single-Level  3  Values  From  HALOE  Monthly  Mean 
Measurements  Over  1996-2005  and  10°S-10°N 


Pressure 

(hPa) 

Case  a: 
Reported 
Errors 

Case  b: 
2*H20/ 
0.5*CH4 

Case  c: 
0.5*H2O/ 
2.0*CH4 

Correlation 

Coefficient 

1.5 

2.87 

2.74 

3.11 

-0.93 

2.2 

2.62 

2.53 

2.73 

-0.96 

3.2 

2.44 

2.37 

2.50 

-0.97 

4.6 

2.57 

2.33 

2.70 

-0.92 
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Table  3.  Single-Level  (3  Values  From  ACE  Occultation 
Measurements  Over  2004-2007  and  10°S-10°N 


Pressure 

(hPa) 

Case  a: 
Reported 
Errors 

Case  b: 
2*H,0/ 
0.5*CH4 

Case  c: 
0.5*H2O/ 
2.0*CH4 

Correlation 

Coefficient 

1.0 

3.54 

2.75 

3.68 

-0.81 

1.5 

2.89 

2.50 

2.99 

-0.89 

2.2 

2.61 

2.23 

2.75 

-0.88 

3.2 

2.30 

2.03 

2.48 

-0.89 

4.6 

2.54 

2.03 

2.92 

-0.81 

absolute  value  of  the  correlation  coefficient  is  at  least  0.65 
in  this  study.  The  (3  values  are  all  statistically  greater  than 
2.0  for  all  pressure  levels  except  for  the  4:1  H20/CH4 
random  error  ratio  case  at  4.6  hPa  and  also  show  an  altitude 
dependence  in  the  same  sense  as  observed  by  HALOE.  (3 
values  were  also  calculated  for  the  overlapping  time  period 
of  HALOE  and  ACE  observations  of  2004-2005.  Good 
agreement  between  the  HALOE  and  ACE  [3  values  were 
found  for  1.5-4. 6  hPa. 

[29]  The  M1PAS  H20  and  CH4  measurements  were 
analyzed  to  calculate  (3  values  for  pressures  between  1 
and  4.6  hPa.  Monthly  mean  residuals  for  H20  and  CH4 
were  calculated  for  the  1 9  months  of  observations  available 
between  September  2002  and  March  2004.  For  the  linear  fit 
calculation,  the  random  uncertainty  for  H20  and  CH4  is 
determined  by  taking  the  root-sum-square  of  the  individual 
random  error  sources  for  the  H20  and  CH4  measurements, 
which  include  detector  noise,  uncertainties  of  parameters  in 
the  retrieval  calculation,  and  line-of-sight  pointing  errors. 
Since  random  errors  are  not  currently  available  for  every 
measurement,  an  individual  profile  of  random  error  sources 
was  used  in  this  analysis.  The  relative  random  errors 
between  the  various  error  sources  for  this  one  profile  should 
be  representative  of  tropical  conditions.  The  calculated  (3 
values  for  each  pressure  level  are  shown  in  Table  4. 
Reasonable  agreement  between  the  MIPAS  (3  values  and 
those  from  HALOE  and  ACE  is  apparent;  however,  the 
monotonic  altitude  dependence  observed  in  the  HALOE 
and  ACE  [3  is  not  seen  in  the  MIPAS  values.  This  discrep¬ 
ancy  could  be  related  to  the  short  time  period  used  for  the 
MIPAS  analysis.  The  actual  random  errors  for  the  1 : 1  H20/ 
CH4  random  error  ratio  case  were  7-10%  for  CH4  and  5— 
10%  for  H20  for  the  range  of  pressures  in  Table  4,  which 
were  reduced  by  the  square  root  of  the  average  number  of 
MIPAS  profiles  per  month  for  the  regression  calculations. 
The  regression  fitting  errors  ranged  from  0.01  to  0.07  for  the 
five  pressures  and  three  random  error  ratio  cases. 

[30]  Having  calculated  a  range  of  f3  values,  we  now 
address  mechanisms  which  may  affect  (3  in  the  tropical 
upper  stratosphere. 

4.1.1.  Propagation  of  Lower/Middle  Stratospheric  HzO 
and  CH4  Variations 

[31]  As  air  from  the  tropical  tropopause  propagates  up¬ 
ward,  we  would  expect  the  tape  recorder  signal  to  decrease 
with  increasing  altitude,  and  indeed  the  H20  variations 
observed  in  the  HALOE  data  generally  decrease  from  68 
to  10  hPa,  and  the  H20  +  2*CH4  variations  decrease  from 
68  to  5  hPa  (see  Figure  3).  But  from  5  hPa  to  1  hPa,  where 
variations  in  CH4  begin  to  become  significant  relative  to 
those  in  H20,  variations  in  H20  and  in  H20  +  2*CH4 
generally  increase.  Under  some  circumstances,  an  increase 


in  temporal  variation  in  H20  +  2*CH4  with  increasing 
altitude  does  occur  in  the  lower  stratosphere,  such  as  the 
increasing  QBO  variations  which  McCormack  and  Siskind 
[2002]  attribute  to  the  interaction  of  QBO-induced  circula¬ 
tion  anomalies  with  the  seasonal  variations  in  H20  mixing 
ratios.  However,  while  some  of  the  variations  in  H20  + 
2*CH4  in  the  upper  stratosphere  are  certainly  driven  by  the 
variations  in  entry  level  H20,  we  argue  that  they  would 
likely  be  small  in  the  tropical  upper  stratosphere  due  to  the 
diminished  amplitude  of  the  tape  recorder  signal  at  1 0  hPa 
compared  to  at  68  hPa  (see  Figure  3a). 

[32]  It  is  not  possible  remove  the  H20  tape  recorder 
signal  from  the  total  H20  variations  in  the  upper  strato¬ 
sphere  due  to  a  blending  of  the  tape  recorder  with  H20 
variations  related  to  seasonal  and  QBO  fluctuations  in 
transport;  however,  we  can  estimate  an  upper  limit  to  the 
magnitude  of  the  H20  tape  recorder  on  f3  in  the  upper 
stratosphere  by  taking  the  annual  amplitude  of  H20  + 
2*CH4  at  10  hPa  and  then  using  this  value  at  pressures 
above  10  hPa.  The  H20  +  2*CH4  signal  is  used  instead  of 
H20  to  estimate  the  H20  tape  recorder  because  a  small 
portion  of  the  annual  cycle  in  H20  at  10  hPa  is  from  an 
annual  transport  cycle,  which  is  canceled  out  in  H20  + 
2*CH4.  This  estimate  will  be  an  upper  limit  because  the 
tape  recorder  signal  should  continue  to  decay  above  1 0  hPa. 
If  we  take  the  tape  recorder  amplitude  at  10  hPa  of  ±0.08 
ppmv  (Figure  3c)  and  subtract  it  from  the  H20  mixing  ratio 
at  1.5,  2.2,  3.2,  and  4.6  hPa,  we  can  estimate  the  maximum 
effect  the  tape  recorder  will  have  on  the  H20/CH4  relation¬ 
ship.  The  H20/CH4  regression  was  carried  out  the  same  as 
before,  but  the  H20  mixing  ratios  were  altered  by  the 
magnitude  of  the  annual  cycle  corresponding  to  the  time 
of  the  year  of  the  measurement.  Since  the  phase  of  the 
annual  cycle  is  also  not  known,  regressions  were  calculated 
for  a  range  of  phases.  The  calculated  f3  value  remained 
above  2.0  for  all  regressions.  For  example,  (3  at  3.2  hPa 
ranged  from  2.51  to  2.24,  and  (3  at  other  pressures  had 
similar  variations  as  a  function  of  the  phase  of  the  annual 
cycle.  Therefore,  we  conclude  that  tape  recorder  effects 
cannot  fully  account  for  the  observed  variations  in  H20  + 
2*CH4  in  the  equatorial  upper  stratosphere,  although  it  is 
possible  that  they  could  slightly  modify  f3  as  a  function  of 
season. 

4.1.2.  Systematic  Data  Uncertainties 

[33]  Systematic  fractional  biases  in  either  H20  or  CH4 
can  change  the  quantity  H20  +  2*CH4  and  its  variations. 
Although  we  have  attempted  to  adjust  the  HALOE  H20  and 
CH4  to  account  for  low  mixing  ratios  relative  to  ACE  and 
other  data  sets,  we  cannot  be  certain  of  the  absolute 
accuracy  of  the  ACE  data  and,  thus,  must  consider  that 


Table  4.  Single-Level  (3  Values  From  MIPAS  Monthly  Mean 
Measurements  Over  September  2002  to  March  2004  and  10°S- 
10°N 


Pressure 

(hPa) 

Case  a: 
Reported 
Errors 

Case  b: 
2*H,0/ 
0.5*01, 

Case  c: 
0.5*H2O/ 
2.0*01, 

Correlation 

Coefficient 

1.0 

3.16 

2.61 

3.29 

-0.86 

1.5 

2.59 

2.26 

2.69 

-0.92 

2.2 

3.27 

3.06 

3.29 

-0.96 

3.2 

2.67 

2.46 

2.71 

-0.94 

4.6 

3.16 

2.79 

3.23 

-0.84 
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Table  5.  Multilevel  f5  Values  From  HALOE  Over  1996-2005, 
ACE  Over  2004-2007,  and  MIPAS  Over  September  2002  to 
March  2004  Between  10°S  and  10°Na 


Pressure 
Range  (hPa) 

HALOE 

ACE 

MIPAS 

1. 5-4.6 

2.65  (-0.93) 

2.48  (-0.95) 

2.77  (-0.98) 

1. 5-2.2 

2.66  (-0.77) 

2.57  (-0.88) 

2.85  (-0.95) 

22-3.2 

2.48  (-0.87) 

2.64  (-0.91) 

2.77  (-0.96) 

3. 2-4. 6 

2.38  (-0.90) 

2.32  (-0.89) 

2.83  (-0.95) 

“Correlation  coefficients  shown  in  parentheses. 

H20  and  CH4  might  still  have  systematic  fractional  uncer¬ 
tainties  which  could  result  in  changes  to  the  variations  in 
H20  +  2*CH4.  In  order  to  reduce  the  variations  of  Et20  + 
2*CEI4,  we  can  infer  from  Figure  2  that  CH4  would  need  to  be 
increased,  EI2O  decreased,  or  some  combination  of  these  two. 

[34]  The  (3  values  of  2.0  could  be  obtained  at  2.2,  3.2,  and 
4.6  hPa  by  either  reducing  H20  by  20%,  increasing  CEI4  by 
25%,  or  a  combination  of  the  two  (e.g.,  —10%  for  EDO, 
+12.5%  for  CH4).  The  intermediate  case  with  ~10% 
systematic  error  for  each  gas  is  comparable  in  value  to  the 
total  reported  systematic  data  uncertainties  [. Harries  et  al., 
1996b;  Park  et  al.,  1996]  hence,  given  that  we  have 
increased  the  ElALOE  H20  mixing  ratios  by  ~5%  to  bring 
them  closer  to  the  ACE  and  ATMOS  values,  this  is  certainly 
not  an  implausible  explanation  for  the  EI2O  +  2*CEl4 
variations  when  we  consider  the  ElALOE  data  only.  How- 
ever,  given  the  preponderance  of  quantitative  evidence 
suggesting  that  the  HALOE  water  vapor  measurements 
are  biased  slightly  low  (see  section  3),  it  seems  unlikely 
that  they  are  actually  too  high.  It  is  also  noted  that  similar 
arguments  can  be  made  regarding  the  effect  of  potential 
systematic  biases  on  8  values  from  the  ACE  and  MIPAS 
data. 

4.1.3.  Vertical  Resolution  Differences 

[35]  Vertical  resolution  differences  between  HALOE  H20 
and  CH4  are  another  factor  that  could  affect  H20  +  2*CH4. 
The  vertical  resolution  of  the  HALOE  H20  measurements  is 
slightly  over  2.0  km  while  the  CH4  measurements  have  a 
4.5  km  resolution.  In  practice,  these  measurements  should 
have  the  same  vertical  resolution  if  they  are  combined  into 
the  quantity  H20  +  2*CH4.  We  note  that  this  issue  does  not 
affect  MIPAS  and  ACE  since  these  H20  and  CH4  profiles 
have  the  similar  vertical  resolutions. 

[36]  We  have  estimated  this  effect  by  smoothing  the  H20 
profiles  to  4.5  km  to  match  the  CH4  vertical  resolution.  This 
smoothing  is  achieved  by  convolving  a  4  km  wide  Gaussian 
function  with  the  H20  profiles  resulting  in  a  net  vertical 
resolution  of  ~4.5  km.  Then,  regression  analyses  were 
carried  out  on  the  individual  H20  and  CH4  measurements 
to  calculate  [3  values.  The  (3  values  were  reduced  from  2.67 
to  2.42  at  1.5  hPa  and  from  2.48  to  2.31  at  2.2  hPa,  for 
example,  but  remained  significantly  above  2.0  throughout 
the  upper  stratosphere. 

4.2.  Multilevel  Correlation  Analysis 

[37]  We  now  consider  the  calculation  of  j3  using  multiple 
pressure  levels  (multilevel  f3  values).  Combining  different 
pressure  levels  has  the  advantage  of  better  covering  the  H20 
and  CH4  variations,  which  will  make  the  (3  calculation  less 
dependent  on  estimated  random  errors.  However,  calculating 
a  (3  by  combining  data  from  different  pressure  levels  does 


make  the  calculation  sensitive  to  any  pressure-dependent 
systematic  error  in  the  measurements.  In  the  absence  of 
measurement  error,  and  provided  the  true  (3  is  similar 
between  levels,  the  two  methods  should  give  very  similar 
results.  For  a  larger  separation  of  the  pressure  surfaces,  the 
H20  and  CH4  changes  represent  the  average  changes  over 
several  pressure  surfaces  due  to  various  transport  processes 
of  possibly  varying  type  and  magnitude,  and  the  calculated 
/ 3  will  have  an  intermediate  value  between  the  single-level  f3 
values  bounding  the  pressure  range. 

[38]  In  Table  5  we  show  (3  values  calculated  between 
several  pressure  surfaces  for  HALOE,  ACE,  and  MIPAS 
data.  For  the  multilevel  analysis,  only  the  data  at  the  two 
pressure  levels  which  define  the  pressure  range  are  used; 
data  for  intermediate  pressure  levels  are  not  included  in  the 
regression.  Figure  4a  shows  the  distribution  of  HALOE 
H20  and  CH4  used  to  calculate  8  between  1.5  and  2.2  hPa 
using  a  4:1  H20/CH4  random  error  ratio.  The  linear  least 
squares  fit  line  to  the  cumulative  distribution  of  data  points 
at  1.5  and  2.2  hPa  is  shown  by  the  solid  black  line.  The 
calculated  f3  of  2.66  between  1.5  and  2.2  hPa  falls  in 
between  the  single-level  (3  values  calculated  using  a  4:1 
H20/CH4  random  error  ratio  at  1.5  and  2.2  hPa,  thus  this 
choice  of  random  errors  and  bias  correction  to  the  HALOE 
data  is  self-consistent.  Figure  4b  shows  the  distribution  of 
ACE  H20  and  CH4  used  to  calculate  /3  between  1.5  and 
2.2  hPa  using  a  1:1  H20/CH4  random  error  ratio.  This 
random  error  ratio  is  used  because  it  leads  to  the  best 
agreement  between  the  ACE  single  and  multilevel  f3  values 
between  1.5  and  4.6  hPa.  A  4:1  H20/CH4  random  error  ratio 
was  used  for  the  MIPAS  (3  value  calculations.  The  ACE  and 
MIPAS  results  are  similar  to  those  from  HALOE.  This 
similarity  is  partially  imposed  because  we  have  bias- 
corrected  the  HALOE  and  MIPAS  profiles  to  match  the 
ACE  data.  The  consistency  between  the  multilevel  and 
single-level  calculations  is,  however,  not  imposed,  hence 
the  general  agreement  between  these  values  and  those  of  the 
single-level  calculations  is  encouraging.  All  the  results 
show  (3  values  above  2.0. 

[39]  In  Figure  5  we  illustrate  the  importance  of  bias 
correcting  the  HALOE  data  to  ACE.  The  linear  least  squares 
fit  lines  shown  here  are  calculated  by  assuming  half  the  CH4 
random  error  and  twice  the  H20  random  error  at  each  level. 
Figure  5a  shows  the  distribution  of  HALOE  H20  versus  CH4 
data  at  1.5  and  4.6  hPa  for  uncorrected  (i.e.,  version  19) 
data.  The  calculated  (3  values  are  3.02  and  2.63  at  1.5  and 
4.6  hPa,  clearly  larger  than  the  value  of  2.10  for  the 
multilevel  case  of  1.5  to  4.6  hPa.  It  is  also  noted  that 
single-level  f3  values  for  intermediate  pressures  between  1.5 
and  4.6  hPa  also  have  values  exceeding  2.0.  Figure  5b 
shows  a  similar  plot  but  for  HALOE  data  which  have  been 
bias  corrected  to  ACE  as  described  in  section  3.  Here,  as 
expected,  the  multilevel  f3  over  the  pressure  range  has  an 
intermediate  value  between  the  two  single-level  f3  values. 
The  similarity  of  the  single  and  multilevel  (3  values  suggests 
that  bias  correcting  the  HALOE  data  to  ACE  has  produced  a 
more  geophysically  consistent  data  set. 

5.  CHEM2D  Modeling 

[40]  Output  from  the  Naval  Research  Laboratory  (NRL) 
zonally  averaged  CHEM2D  model  has  been  used  to  inves- 
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Figure  5.  Water  vapor  versus  methane  at  1.5  and  4.6  hPa  measured  by  HALOE  from  1996  to  2005 
using  (a)  version  19  data  and  (b)  data  bias  corrected  to  ACE.  The  data  points  and  linear  least  squares  fit 
line  to  the  data  are  shown  in  blue  for  1.5  hPa  and  in  red  for  4.6  hPa.  The  linear  least  squares  fit  line 
through  the  cumulative  distribution  of  points  at  1.5  and  4.6  hPa  is  shown  in  black.  All  data  shown  are  for 
the  latitude  range  of  10°S  -  10°N.  The  slopes  of  the  linear  least  squares  fit  lines  are  shown  in  the  top  right- 
hand  comer  of  each  plot.  The  dashed  line  has  a  slope  of  —2.0. 


tigate  the  hydrogen  chemistry  of  the  upper  stratosphere. 
CHEM2D  is  a  global  middle  atmosphere  model  extending 
from  the  surface  to  ~120  km  altitude  with  interactive 
radiative,  photochemical  and  dynamical  schemes.  The  ver¬ 
sion  of  CHEM2D  used  in  the  present  study  has  88  vertical 
levels  with  a  vertical  grid  spacing  of  1.5  km  [ McCormack  et 
al.,  2007].  CHEM2D  has  been  shown  to  successfully 
simulate  the  general  distribution  and  variability  of  tracers 
such  as  H20,  CH4,  and  N20  [ Summers  et  al.,  1997; 
McCormack  and  Siskind,  2002].  Here  we  analyze  monthly 
zonal  mean  H20,  CH4  and  H2  fields  from  a  30  year 
CHEM2D  simulation  in  which  the  entry  value  of  H20 
mixing  ratio  at  the  tropical  tropopause  (17  km)  varies 
sinusoidally  with  an  amplitude  of  1.5  ppmv  about  a  mean 
value  of  3.5  ppmv.  No  trends  or  other  interannual  variations 
in  H20  entering  the  middle  atmosphere  from  the  tropical 
troposphere  are  present.  Model  CH4  mixing  ratios  are  held 
fixed  at  1.75  ppmv  in  the  troposphere  throughout  the 
simulation.  The  model  includes  an  interactive  parameteri¬ 
zation  of  the  QBO  in  equatorial  stratospheric  winds.  The 
QBO  modulates  the  strength  of  the  modeled  tropical  up- 
welling  in  the  lower  stratosphere  and  produces  a  small  QBO 
in  model  H20  and  CH4  [ McCormack  and  Siskind,  2002], 
[4i]  Figure  6  shows  H20,  CH4,  H20  +  2*CH4,  H2,  and 
H20  +  2*CH4  +  H2  at  a  pressure  of  1 .4  hPa  for  years  20-25 
of  the  simulation  averaged  over  10°S  to  10°N.  While  the 
variability  in  H20  and  CH4  is  somewhat  different  than  that 
seen  in  the  HALOE  data,  the  CHEM2D  model  shows  that 
H20  +  2*CH4  has  time  variations  which  are  anticorrelated 
with  the  CH4  variations  in  the  same  sense  as  the  observa¬ 
tions.  As  expected,  the  H20  +  2*CH4  variations  are  also 
anticorrelated  with  variations  in  H2,  such  that  H20  + 
2*CH4  +  H2  is  nearly  constant  in  time.  The  remaining 
fluctuations  in  H20  +  2*CH4  +  H2  are  due  to  seasonal 
variations  in  the  H20  tape  recorder  signal;  therefore,  we 
have  used  this  time  series  to  estimate  the  tape  recorder 
signal  by  calculating  the  residuals  from  the  mean  H20  + 
2*CH4  +  H2  mixing  ratio.  This  tape  recorder  estimate  is 
subtracted  from  the  H20  time  series  at  each  pressure,  and 


regression  fits  are  made  to  the  H20  versus  CH4  changes. 
The  calculated  single-level  p  values  based  on  the  H20  and 
CH4  variations  at  1.2,  1.4,  2.1,  3.0,  and  4.4  hPa  for  years 
20-25  are  2.42,  2.38,  2.29,  2.20,  and  2.1 1,  respectively.  The 
calculated  multilevel  /3  value  between  1.2  and  4.4  hPa  is 
2.27,  or  roughly  an  average  of  the  single-level  P  values  over 
this  pressure  range.  These  values  are  all  greater  than  2,  but 
they  are  somewhat  smaller  than  the  values  determined  from 
observations.  We  would  expect  that  the  P  values  are  a 
function  of  both  circulation  and  chemistry,  hence  even  if 
the  model  chemistry  is  correct,  it  is  possible  that  the 
disagreement  in  the  model  and  observed  (3  values  may  be 
attributed  to  differences  between  the  CHEM2D  circulation 
and  the  actual  circulation.  The  single-level  P  values  also 
show  a  clear  altitude  dependence  after  the  effects  of  the  H20 
tape  recorder  were  removed  in  a  similar  sense  to  behavior  of 
the  HALOE  and  ACE  p  values.  A  greater  proportion  of  H2 
loss  versus  production  is  likely  responsible  for  most  of  the 
altitude  dependence  though  some  may  also  be  due  other 
transport  processes  not  properly  simulated  in  the  model. 


Year 

Figure  6.  Time  series  of  water  vapor,  methane,  molecular 
hydrogen,  H20  +  2*CH4,  and  H20  +  2*CH4  +  H2  at  1.4  hPa 
averaged  over  10°S-10°N  from  CHEM2D.  The  time  series 
shown  is  for  years  20-25  of  a  30  year  simulation.  Note  that 
some  of  the  time  series  have  been  shifted  to  fit  the  plot. 
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Figure  7.  Mean  profiles  of  water  vapor,  methane,  molecular  hydrogen,  H20  +  2*CH4,  and  H20  + 
2*CH4  +  H2  from  (a)  CHEM2D  averaged  over  years  20-25  of  a  30  year  simulation,  (b)  HALOE 
averaged  over  1996-2005,  (c)  ACE  averaged  over  2004-2007,  and  (d)  MIPAS  averaged  over  October 
2002  to  September  2003.  All  profiles  were  averaged  over  10°S—  10°N,  and  the  water  vapor,  methane,  and 
molecular  hydrogen  profiles  were  shifted  to  fit  the  plots.  The  line  description  for  each  profile  is  shown  in 
Figure  7a.  We  note  that  CHEM2D  simulates  each  component  of  H20  +  2*CH4  +  H2  explicitly,  while  H2 
is  inferred  from  the  observations  using  equation  (3). 


[42]  Mean  profiles  of  the  relevant  hydrogen  parameters 
from  CHEM2D  averaged  over  years  20-25  and  10°S  to 
10°N  are  shown  in  Figure  7a.  Similarly  calculated  profiles 
from  the  HALOE,  ACE,  and  MIPAS  observations  shown  in 
Figures  7b-7d  will  be  discussed  later.  As  expected  from  the 
CHEM2D  single-level  (3  values  >2,  a  slight  positive  vertical 
gradient  is  observed  in  H20  +  2*CH4  above  ^5  hPa.  The 
increase  in  model  H20  +  2*CH4  of  about  0.15  ppmv 
between  the  5  and  1.5  hPa  levels  is  balanced  by  a  compa¬ 
rable  decrease  in  H2,  suggesting  that  at  least  0.15  ppmv 
more  H2  is  lost  through  oxidation  than  produced  from 
CH20  photolysis  in  the  equatorial  upper  stratosphere. 

6.  Measurement  Time  Series  and  Vertical  Profiles 

[43]  The  above  results  highlight  the  nonconservation  of 
H20  +  2*CH4  in  the  equatorial  upper  stratosphere  which 
results  from  H2  not  being  in  steady  state.  As  expected, 
model  results  from  CHEM2D  show  H20  +  2*CH4  varia¬ 
tions  which  are  balanced  primarily  by  H2  variations,  so  that 
total  hydrogen,  Hy,  is  conserved  by  adding  H2  to  H20  + 
2*CH4, 

Hy  =  H20  +  2*CH4  +  H2 .  (2) 


Variations  in  total  hydrogen  can  then  be  calculated  as  a 
function  of  pressure  and  time  from  the  observational  data  by 
relating  the  H2  term  in  equation  (2)  with  the  variations  in 
H20  +  2*CH4.  We  believe  that  this  is  a  valid  approach  in 
light  of  the  CHEM2D  results  and  having  ruled  out  other 
potential  sources  of  the  H20  +  2*CH4  variations.  If  the  H2 
term  is  expanded  into  a  mean  and  residual  term  and 
equation  (1)  is  substituted  in  (2),  then  Hy  at  pressure  level  p; 
can  be  expressed  solely  in  terms  of  H20,  CH4,  and  the  3 
values  (along  with  an  assumed  constant  H2  mixing  ratio  of 
0.55  ppmv  at  5  hPa  (level  1  in  (3)), 

Hy(Pi,  t)  =  H20(Pi,t)  +  2*CH4(Pi,t)  +  <H2(7  hPa)} 

+  £  {[0.5*(/3(Pi)  +  fe,))  -  2]*[(CH4(Pi)) 
i=ij 

-  (CH4(pi_1))]}  +  (/3(Pi)  -  2)*dCELt(Pi,  t).  (3) 

The  third  term  of  equation  (3)  is  the  long-term  average  H2 
mixing  ratio  at  5  hPa  (taken  from  in  situ  measurements 
showing  average  H2  mixing  ratios  of  0.55  ppmv  in  the  8— 
10  hPa  range  [Zoger  et  al.,  1999;  Rohs  et  al.,  2006]),  while 
the  fourth  term  is  the  summation  of  the  changes  in  the  long¬ 
term  average  CH4  mixing  ratio  between  5  hPa  and  the 
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Figure  8.  Time  series  of  corrected  HALOE  (black  line),  ACE  (blue  line),  and  corrected  MIPAS  (red 
line)  (top)  water  vapor,  two  times  methane,  and  inferred  molecular  hydrogen  and  (bottom)  H20  +  2*CH4 
and  Hy  at  3.2  hPa  averaged  over  10°S-10°N.  The  H20  and  H2  time  series  were  shifted  to  fit  the  plot. 
Data  markers  on  the  top  of  each  plot  show  the  times  of  monthly  observations  (at  least  12  profiles  required 
per  month)  for  each  instrument  (same  color  connotations). 


pressure  of  interest  multiplied  times  (/3  —  2)  where  f3  is  an 
average  of  the  value  at  the  chosen  pressure  and  one  level 
below.  We  choose  5  hPa  as  the  “reference”  pressure  in  the 
calculation  since  H20  +  2*CH4  has  minimal  variations  near 
this  pressure  and  because  this  was  the  highest  pressure  at 
which  p  calculations  were  possible.  For  the  fifth  tenn,  we 
have  calculated  CH4  monthly  residuals  from  the  overall 
time  mean  for  a  particular  instrument.  Finally,  the  p  values 
in  the  fourth  term  are  pressure  dependent  between  5  hPa  and 
the  chosen  pressure,  while  the  P  value  in  the  fifth 
corresponds  to  the  pressure  at  which  Hy  is  calculated. 

[44]  Although  time  series  of  Hy  can  be  calculated  using 
the  p  values  calculated  from  the  correlation  analysis,  in  this 
calculation,  we  search  for  a  value  of  /3  such  that  Hy  at  each 
level  is  conserved.  The  FIALOE  data  were  used  to  calculate 
the  p  value  that  minimizes  the  variations  in  Hy  (hereafter 
termed  the  “residual”  p  value)  since  these  data  have  the 
longest  record  of  simultaneous  H20  and  CFl4  observations 
in  this  region  of  the  tropical  upper  stratosphere  and,  thus, 
should  provide  the  best  long-term  average  3  value.  The  time 
series  were  again  restricted  to  1996-2005  in  order  to  avoid 
a  possible  influence  of  the  early  1990s  increase  in  upper 
stratospheric  FI20  on  the  residual  3  value  calculations. 
Multiple  time  series  of  FIy  were  calculated  at  each  pressure 
level  (1.5,  2.2,  3.2,  and  4.6  hPa)  using  equation  (3)  where 
the  p  value  was  varied  from  2.0  to  4.0  in  increments  of  0.02. 
Then,  standard  deviations  of  the  Hy  time  series  at  each 
pressure  level  were  calculated  in  order  to  find  the  p  (at  each 
pressure)  that  gave  the  smallest  standard  deviation.  The 
calculated  residual  P  values  were  2.66,  2.56,  2.38,  and  2.30 
for  1.5,  2.2,  3.2,  and  4.6  hPa,  respectively.  These  P  values 
result  in  a  reduction  of  the  standard  deviation  in  FIy  of  4%  at 
4.6  hPa  up  to  a  maximum  of  19%  at  2.2  hPa  compared  to 


time  series  of  Hy  computed  with  a  3  value  of  2.  Time  series 
of  FIy  for  the  FIALOE,  ACE,  and  MIPAS  data  are  then 
calculated  at  the  four  pressure  levels  using  the  HALOE 
derived  residual  P  values. 

[45]  Figures  8  shows  time  series  of  Hy  at  3.2  hPa 
averaged  over  10°S  to  10°N  along  with  time  series  of 
H20,  2*CH4,  inferred  H2,  and  H20  +  2*CH4  where  the 
H2  mixing  ratios  are  calculated  from  the  last  three  terms  of 
equation  (3).  The  Hy  and  H20  +  2*CH4  time  series  are 
presented  with  a  smaller  vertical  scale  in  order  to  better 
show  the  small  variations.  Figure  8  shows  that  the  variations 
seen  in  Hy  are  noticeably  reduced  compared  with  those  in 
the  H20  +  2*CH4  time  series.  The  ACE  and  MIPAS  Hy 
time  series  also  show  an  overall  reduction  in  total  variability 
when  the  residual  3  values  HALOE  are  used  and  overlap 
with  the  HALOE  Hy  time  series.  The  overall  agreement  in 
Hy  between  the  three  instruments  is  very  encouraging  and 
provides  confidence  in  the  mixing  ratio  of  total  hydrogen 
and,  as  shown  later,  stratospheric  entry  level  H20.  These 
results  highlight  the  importance  of  understanding  variations 
in  H2  as  well  as  those  of  H20  and  CH4  in  this  region  of  the 
atmosphere. 

[46]  Least  squares  regression  fits  were  done  to  the 
HALOE  Hy  time  series  over  1996-2005  averaged  over 
10°S  to  10°N  as  was  done  in  Figure  3  for  H20,  CH4,  and 
H20  +  2*CH4.  Figure  9  shows  the  amplitudes  of  the  QBO, 
annual,  and  semiannual  cycles  as  a  function  of  pressure 
calculated  from  the  fit.  Since  Hy  is  only  calculated  for 
pressures  less  than  5  hPa,  the  amplitudes  for  H20  +  2*CH4 
are  shown  below  this  pressure  level.  A  comparison  of  the 
cycle  amplitudes  of  H20  +  2*CH4  (Figure  3c)  to  those  of 
Hy  (Figure  9)  reveals  that  the  increases  in  the  cycle 
amplitudes  above  5  hPa  seen  in  H20  +  2*CH4  are  consid- 
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Figure  9.  Profiles  of  annual  (solid  black  line),  semiannual 
(dashed  red  line),  and  QBO  (dash-dotted  blue  line)  cycle 
amplitudes  from  a  least  squares  fit  to  HALOE  Hy  time  series 
from  1996  to  2005  averaged  over  10°S-10°N.  Below  5  hPa 
the  cycle  amplitudes  are  calculated  from  the  H20  +  2*CH4 
time  series  instead  of  Hy. 


erably  reduced  in  Hy.  For  instance,  at  2.2  hPa,  the  QBO 
amplitude  in  Hy  was  reduced  by  66%  and  the  annual  cycle 
amplitude  was  reduced  by  52%  relative  to  the  respective 
amplitudes  in  the  H20  +  2*CH4  time  series.  A  general 
decrease  in  the  cycle  amplitudes  with  decreasing  pressure  is 
now  observed  over  the  entire  stratosphere,  with  the  excep¬ 
tion  near  15  hPa  where  the  QBO  winds  and  H20  tape 
recorder  signals  interact  to  produce  a  local  maximum  in  the 
QBO  amplitude. 

[47]  If  the  time  series  of  Hy  has  successfully  removed  the 
variations  in  H20  +  2*CH4  which  are  related  to  upper 
stratospheric  chemistry,  then  the  remaining  variations  result 
from  variations  in  H20  and  CH4  entering  the  stratosphere  in 
the  tropics.  We  can  therefore  use  this  time  series  along  with 
CH4  time  series  and  the  residual  (3  values  to  calculate  the 
time  series  of  stratospheric  entry  H20  mixing  ratios.  In 
order  to  obtain  the  entry  H20  from  Hy,  we  subtract  the  H2 
mixing  ratio  and  twice  the  CH4  mixing  ratio  at  98  hPa.  For 
CH4  at  the  tropopause,  the  2*CFl4  mixing  ratio  between 
1999  and  2005  (avoiding  the  immediate  post-Pinatubo 
period)  is  3.58  ppmv.  For  H2  at  the  tropopause,  a  value  of 
0.55  ppmv  is  used  [Rohs  et  al.,  2006].  If  we  use  the  3.2  hPa 
Hy  time  series,  we  can  derive  the  stratospheric  entry  H20 
result  shown  in  Figure  10.  Because  it  is  based  upon  the 
3.2  hPa  time  series,  the  calculated  entry  H20  time  series  at 
98  hPa  is  of  course  both  smoothed  and  delayed  by  several 
years  compared  to  the  actual  entry  level  time  series.  It 
should,  however,  be  suitable  for  those  interested  in  the 
average  H20  entry  level  values  averaged  over  multiyear 
periods.  For  example,  comparisons  between  the  calculated 
entry  H20  time  series  and  long-term  records  of  balloon- 
borne  midlatitude  H20  measurements  in  the  lower  strato¬ 
sphere  [Oltmans  et  al.,  2000;  Scherer  et  al.,  2008]  would  be 
appropriate. 

[48]  The  uncertainty  in  the  stratospheric  entry  H20  mixing 
ratio  is  driven  by  the  errors  in  H20,  2*CH4,  and  (8  —  2)*d 
(CH4)/dp  in  the  Hy  calculation.  Since  the  HALOE  data  have 
been  bias  corrected  to  the  ACE  data,  the  H20  and  CH4 


systematic  errors  are  effectively  equal  to  those  for  ACE 
H20  (5%)  and  CH4  (10%)  in  the  upper  stratosphere.  This 
5%  uncertainty  in  H20  at  3.2  hPa  introduces  —0.28  ppmv 
uncertainty  in  entry-level  H20,  while  the  10%  uncertainty  in 
CH4  introduces  an  uncertainty  of  — 0. 1 8  ppmv.  The  error  in 
8,  estimated  from  the  difference  in  the  HALOE  single-level 
8  values  between  the  1:1  and  4:1  H20/CH4  random  error 
ratio  cases  for  the  monthly  averaged  data  (see  Table  2),  is 
less  than  0. 1  at  most  pressures.  The  contribution  of  uncer¬ 
tainties  in  (3  to  the  total  error  in  the  entry  H20  mixing  ratio 
is  therefore  —0.02  ppmv.  The  total  uncertainty  in  the  entry 
H20  mixing  ratio  from  the  above  three  error  sources  is 
calculated  to  be  —0.3  ppmv. 

[49]  Using  the  HALOE  derived  residual  /3  values  which 
minimized  Hy  variations  at  each  level,  we  then  investigated 
the  vertical  profile  distributions  of  CH4,  H20  +  2*CH4, 
inferred  H2,  and  Hy  in  the  equatorial  upper  stratosphere 
from  HALOE,  ACE,  and  MIPAS.  Multiyear  profiles  of 
CH4,  H20  +  2*CH4,  H2,  and  Hy  averaged  over  10°S  to 
10°N  for  HALOE,  ACE,  and  MIPAS  data  are  shown  in 
Figures  7b-7d.  The  mean  profiles  for  HALOE  and  ACE  are 
calculated  by  first  taking  multiyear  monthly  averages,  then 
averaging  the  monthly  averages.  For  MIPAS,  a  yearly 
average  profile  was  calculated  from  monthly  zonal  mean 
data  from  October  2002  to  September  2003.  Results  from 
MIPAS  are  similar  using  other  1  year  periods.  Yearly 
averages  were  calculated  with  all  three  data  sets  so  as  to 
remove  potential  seasonal  cycle  H20  and  CH4  variations. 
The  solid  H20  +  2*CH4  profiles  between  1.5  and  5  hPa  do 
not  show  the  typical  near  zero  vertical  gradient  seen  in  the 
extratropical  lower  stratosphere  [Abbas  et  al.,  1996;  Engel 
et  al.,  1996;  Hurst  et  al.,  1999;  Nassar  et  al,  2005],  but 
instead  a  slight  positive  value  with  mixing  ratios  increasing 
by  about  0.4  ppmv  over  this  range.  We  note  that  a  previous 
study  of  MIPAS  data  [Juckes,  2007]  has  revealed  this 
feature  in  H20  +  2*CH4.  The  inferred  profiles  of  H2  for 
each  instrument  are  also  shown  (the  third  and  fourth  terms 
of  equation  (3)),  which  shows  a  —0.3  ppmv  reduction  above 
5  hPa.  Mean  profiles  of  Hy  are  calculated  by  summing  the 
mean  H20  +  2*CH4  and  H2  profiles.  We  note  that  the  fifth 
term  of  equation  (3)  does  not  enter  the  calculation  for  Hy 
because  the  CH4  residuals  average  to  zero  in  the  time  mean. 
The  Hy  profiles  generally  show  that  the  8  values  which 
minimize  the  time  variation  in  Hy  on  individual  levels  also 


Figure  10.  Time  series  of  HALOE  stratospheric  entry 
water  vapor  at  98  hPa  averaged  over  10°S-10°N  and 
derived  from  the  HALOE  Hy  time  series  at  3.2  hPa.  Note 
that  this  time  series  is  both  highly  smoothed  and  delayed 
several  years  from  the  actual  entry  water  vapor  time  series. 
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Figure  11.  Vertical  profiles  of  the  change  in  H20,  H20  + 
2*CH4,  and  Hy  per  year  for  1992-1996  averaged  over 
10°S-10°N.  The  yearly  changes  were  calculated  from  a 
linear  least  squares  fit  to  the  HALOE  time  series  which 
included  QBO,  annual,  semiannual,  and  a  linear  trend  term. 

reduce  the  variation  of  total  hydrogen  with  altitude.  This  is 
perhaps  the  most  compelling  evidence  that  the  systematic 
adjustment  which  we  have  made  to  the  HALOE  data  based 
on  comparisons  with  ACE  is  reasonable. 

[50]  The  general  altitude  dependence  of  the  mean  H20  + 
2*CH4  and  H2  profiles  is  similar  between  the  CHEM2D 
results  (Figure  7a)  and  the  observations  (Figures  7b-7d), 
though  the  vertical  gradients  are  about  half  as  large  for  the 
model.  As  mentioned  earlier,  we  suspect  that  discrepancies 
between  the  modeled  and  actual  atmospheric  circulation 
could  contribute  to  these  differences.  Indeed,  a  slight 
difference  in  circulation  between  the  model  and  observa¬ 
tions  is  supported  by  the  slightly  larger  drop  in  CH4  with 
altitude  in  the  CHEM2D  model  relative  to  the  CH4  meas¬ 
urements  shown  in  Figure  7. 

[51]  An  interesting  ramification  of  the  Hy  calculation  is 
that  the  increase  observed  in  H20  +  2*CH4  during  the  early 
and  middle  1990s  [e.g.,  Nedoluha  et  al.,  2003]  is  slightly 
reduced  in  magnitude.  We  analyze  the  yearly  changes  by 
performing  least  squares  fits  to  the  H20,  H20  +  2*CH4,  and 
Hy  time  series  over  the  1992-1996  time  period  where  we 
have  included  a  linear  trend  term  to  the  fit  in  addition  to 
QBO,  annual,  and  semiannual  terms.  The  yearly  changes  for 
each  parameter  are  derived  from  the  linear  fit  term,  which  is 
plotted  at  4  pressure  levels  (Figure  11).  This  shows  that  at 
3.2  hPa  about  40%  of  the  total  change  in  H20  is  attributable 
to  changes  in  CH4,  and  ~  1 2%  of  the  remaining  variation  is 
attributable  to  changes  in  H2.  The  1  sigma  uncertainties  as 
calculated  from  the  trend  tenn  in  the  fits  are,  at  2.2  hPa, 
~0.03  ppmv/yr  for  H20,  ~0.02  ppmv/yr  for  CH4,  and 
~0.01  ppmv/yr  for  Hy.  The  trends  are  highly  correlated 
over  these  levels,  so,  while  these  uncertainties  are  larger 
than  the  level-to-level  differences,  the  altitude  dependence 
may  nevertheless  be  significant.  Figure  11  reveals  an 
altitude  dependence  in  the  yearly  changes  for  H20.  Most, 
but  not  all,  of  this  altitude  dependence  is  removed  in  H20  + 
2*CH4,  and  the  altitude  dependence  is  further  reduced  in  the 
Hy  profile.  The  Hy  change  profile  is  physically  reasonable 
since  it  implies  that  the  change  in  entry  level  H20  is 
affecting  all  levels  in  the  equatorial  upper  stratosphere 
nearly  equally,  but  with  a  very  slightly  smaller  increase  at 
the  highest  level.  The  similarity  of  Hy  changes  throughout 
the  upper  stratosphere  is  plausible  given  the  short  time  scale 
of  circulation  in  the  upper  stratosphere  of  ~  1  year  [Hall  and 


Waugh,  1997]  (compared  with  the  1992-1996  time  period) 
and  the  near  linearity  of  the  early  and  middle  1990s  Hy 
“trends”  (see  Figure  8). 

[52]  The  early  to  middle  2000s  time  period  was  examined 
to  investigate  possible  changes  in  H20  in  the  equatorial 
upper  stratosphere  given  the  documented  drying  of  the 
equatorial  lower  stratosphere  in  2001  [Randel  et  al., 
2006],  A  similar  analysis  to  that  performed  for  the  early 
1990s  was  done  for  this  time  period,  but  the  upper  strato¬ 
spheric  changes  in  H20,  H20  +  2*CH4,  and  Hy  were  much 
smaller  than  in  the  early  1990s  and  no  clear  conclusions 
could  be  drawn  for  this  period. 


7.  Discussion  and  Conclusions 

[53]  We  have  presented  data  from  the  HALOE,  ACE,  and 
MIPAS  instruments  which  reflect  the  nonconservation  of 
H20  +  2*CH4  on  individual  pressure  surfaces  over  time  or 
as  a  time-averaged  quantity  as  a  function  of  pressure  in  the 
equatorial  upper  stratosphere.  Several  plausible  mechanisms 
were  investigated  to  explain  the  variations  in  H20  +  2*CH4 
such  as  systematic  data  biases  and  vertical  resolution  differ¬ 
ences  in  the  HALOE  H20  and  CH4  data.  The  impacts  of 
these  mechanisms  were  found  to  be  too  small  to  account  for 
the  size  of  the  H20  +  2*CH4  changes. 

[54]  The  equatorial  upper  stratosphere  is  the  ideal  region 
to  calculate  f3  values  since  the  CH4  changes  are  large,  the 
descent  of  air  from  higher  altitudes  is  small,  and  influence 
of  variations  in  stratospheric  entry  level  H20  is  much 
smaller  than  in  the  lower  stratosphere.  After  applying  bias 
corrections  to  HALOE  and  MIPAS  H20  and  CH4  profiles,  13 
values  calculated  from  the  HALOE,  ACE,  and  MIPAS 
showed  good  agreement  both  on  single  pressure  surfaces 
and  over  a  range  of  pressures  in  the  equatorial  upper 
stratosphere  with  values  considerably  above  2.  We  calcu¬ 
lated  values  of  f3  which  minimized  the  variability  of  Hy,  or 
total  hydrogen,  time  series  on  individual  pressure  surfaces. 
The  smoother  Hy  time  series  (compared  with  the  H20  + 
2*CH4  time  series)  was  then  used  to  derive  a  time  series  of 
stratospheric  entry  level  H20  near  the  tropopause.  The 
calculation  of  the  stratospheric  entry  H20  mixing  ratio 
based  on  Hy  in  the  tropical  upper  stratosphere  provides  a 
more  appropriate  value  compared  to  one  based  on  H20  + 
2*CH4  due  to  the  necessity  in  accounting  for  the  net 
reduction  in  H2  between  the  tropopause  and  the  upper 
stratosphere  at  tropical  latitudes. 

[55]  The  measurements  are  supported  by  atmospheric 
simulations  by  the  CHEM2D  model  which  show  not  only 
variations  in  H20  +  2*CH4  in  the  upper  stratosphere  but 
also  variations  in  H2  of  similar  magnitude  but  anticorrelated 
in  time  so  that  H20  +  2*CH4  +  H2  is  nearly  conserved. 
Theoretical  analysis  by  Le  Texier  et  al.  [1988]  has  shown 
that  oxidation  of  H2  plays  a  key  role  in  the  H20  budget  of 
the  stratosphere.  Le  Texier  et  al.  [1988]  showed  that, 
between  40  and  50  km,  there  was  net  loss  of  H2  which 
drives  additional  H20  production  thus  producing  positive 
vertical  gradients  in  H20  +  2*CH4.  This  behavior  is  what  is 
seen  in  the  HALOE,  ACE,  and  MIPAS  profiles  of  H20  + 
2*CH4  and  inferred  H2  in  the  upper  stratosphere,  where 
CH4  oxidation  occurs  most  efficiently.  The  amount  of  the 
excess  H20  +  2*CH4  above  5  hPa  is  ~0.4  ppmv,  which  is 
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nearly  offset  by  a  ~0.3  ppmv  net  reduction  of  H2  (see 
Figure  7). 
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